On a complex differential Riccati equation 
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Abstract 

We consider a nonlinear partial differential equation for complex- 
, valued functions which is related to the two-dimensional stationary 

Schrodinger equation and enjoys many properties similar to those of 
the ordinary differential Riccati equation such as the famous Euler the- 
[-^. . orems, the Picard theorem and others. Besides these generalizations 

of the classical "one-dimensional" results we discuss new features of 
the considered equation including an analogue of the Cauchy integral 
. theorem. 

o 



1 Introduction 



The ordinary differential Riccati equation 

u' = pu 2 + qu + r (1) 

has received a great deal of attention since a particular version was first 
studied by Count Riccati in 1724, owing to both its specific properties and 
the wide range of applications in which it appears. For a survey of the 
history and classical results on this equation, see for example [7] and [25J. 
This equation can be always reduced to its canonical form (see, e.g., [6], [10]). 

y' + y 2 = v, (2) 
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and this is the form that we will consider. 

One of the reasons for which the Riccati equation has so many applica- 
tions is that it is related to the general second order homogeneous differential 
equation. In particular, the one-dimensional Schrodinger equation 



This substitution, which as its most spectacular application reduces Burger's 
equation to the standard one-dimensional heat equation, is the basis of the 
well-developed theory of logarithmic derivatives for the integration of non- 
linear differential equations [19]. A generalization of this substitution will 
be used in this work. 

A second relation between the one-dimensional Schrodinger equation and 
the Riccati equation is as follows. The one- dimensional Schrodinger operator 
can be factorized in the form 

if and only if (j2J) holds. This observation is a key to a vast area of research 
related to the factorization method (see, e.g., [22J and [23]) and to the Dar- 
boux transformation (see, e.g., [20], [2j] and [26]). In the present work we 
consider a result similar to ([3]) but already in a two-dimensional setting. 

Among the peculiar properties of the Riccati equation stand out two 
theorems of Euler, dating from 1760. The first of these states that if a 
particular solution y of the Riccati equation is known, the substitution y = 
y + z reduces (j2J to a Bernoulli equation which in turn is reduced by the 
substitution z = - to a first order linear equation. Thus given a particular 
solution of the Riccati equation, it can be linearized and the general solution 
can be found in two integrations. The second of these theorems states that 
given two particular solutions yo, y\ of the Riccati equation, the general 
solution can be found in the form 



d 2 u + vu = 



(3) 



is related to (j2J by the easily inverted substitution 



y = — 



y = 



ky exp(Jy - y{) - y 1 
kexp(Jy -y x )-l 



(5) 
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where k is a constant. That is, given two particular solutions of d2J), the 
general solution can be found in one integration. 

Other interesting properties are those discovered by Weyr and Picard 
([29J, [7J). The first is that given a third particular solution 2/3, the general 
solution can be found without integrating. That is, an explicit combination 
of three particular solutions gives the general solution. The second is that 
given a fourth particular solution y 4 , the cross ratio 

(j/i -3/2X2/3-2/4) 

(yi -2/4) (j/3 -2/2) 

is a constant. Thus the derivative of this ratio is zero, which holds if and 
only if the numerator of the derivative is zero, that is, if and only if 

(yi-y4)(y3-y2)((yi-y2)(y3-y4))'-(yi-y2)(y3-y4)((yi-y4)(y3-y2))' = 0. 

Dividing by {y t - 2/2X2/3 ~ V^ivi ~ 2/4) (2/3 - 2/2), we see that Picard's 
theorem is equivalent to the statement 

(2/1 ~ 2/2)' (2/3 ~ 2/4) 7 _ (2/1 - 2/4)' _ (2/3 ~ 2/2)' = Q , g v 
2/i - 2/2 2/3 - 2/4 2/i - 2/4 2/3-2/2 

In the present work we study the following equation 

dzQ + \Q\ 2 = v (7) 

where z is a complex variable, = \(d x + id y ), Q is a complex valued 
function of z and f is a real valued function. Note that this equation is dif- 
ferent from the complex Riccati equation studied in dozens of works where 
it is supposed to have the form ([T]) with complex analytic coefficients p, q 
and r and a complex analytic solution u (see, e.g., [9]). We do not sup- 
pose analyticity of the complex functions involved and show that equation 
((7j) unlike the equation considered in [9] is related to the two-dimensional 
stationary Schrodinger equation in a similar way as the ordinary Riccati and 
Schrodinger equations are related in dimension one. Moreover, we establish 
generalizations of the Euler and Picard theorems and obtain some other re- 
sults which are essentially two-dimensional, for example, an analogue of the 
Cauchy integral theorem for solutions of the complex Riccati equation (|7j). 

Equation © first appeared in [H] as a reduction to a two-dimensional 
case of the spatial factorization of the stationary Schrodinger operator which 
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was studied in a quaternionic setting in [TJ, [3J, [H], [T2], [13], P3] and later 
on using Clifford analysis in [2] and [8]. 

The ordinary Riccati equation is at the heart of many analytical and nu- 
merical approaches to problems involving the one-dimensional Schrodinger 
and Sturm-Liouville equations. Here we furnish a complete structural anal- 
ogy between dimensions one and two regarding the relationship between the 
Schrodinger and the Riccati equations. Besides, the deep similarity between 
the ordinary Riccati equation and ([7]) strongly suggests that many known 
applications of the ordinary Riccati equation can be generalized to the two- 
dimensional situation and many new aspects such as theorem [12] will become 
manifest . 



2 Some preliminary results on the stationary 
Schrodinger equation and a class of pseu- 
doanalytic functions 

We need the following definition. Consider the equation d-^p = $ in the 
whole complex plane or in a convex domain, where $ = $x + i$ 2 is a given 
complex valued function whose real part $i and imaginary part $2 satisfy 
the equation 

- d x § 2 = 0, (8) 

then we can reconstruct ip up to an arbitrary real constant c in the following 

way 

ip(x,y) = 2[f $ 1 (r},y)dr J + j $ 2 (x , £)d£ ) + c (9) 

\Jxo Jyo J 

where (xo,yo) is an arbitrary fixed point in the domain of interest. 
By A we denote the integral operator in ([9]): 

A[$](x,y) = 2 ( j <S> 1 ( v ,y)dri+ [ $ 2 (x , c. 

\Jxo Jyo J 

Note that formula ([9]) can be easily extended to any simply connected domain 
by considering the integral along an arbitrary rectifiable curve Y leading from 

(xo,yo) to (x,y), 

<p(x, y) = 2 ^ $icfe + $ 2 dy J + c. 
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Thus if $ satisfies (jSJ), there exists a family of real valued functions ip such 
that dzf = given by the formula <p = A[<f>]. In a similar way we introduce 
the operator 

y) = 2 (^j <&idx - <& 2 dy*j + c 
which is applicable to complex functions satisfying the condition 



d y Q x + <9 X .$ 2 = (10) 

and corresponds to the operator d z . 

Let / denote a positive twice continuously differentiable function defined 
in a domain ficC. Consider the following Vekua equation 

Wz = jW in ft (11) 

where the subindex z means the application of the operator (analogously, 
with the aid of the subindex z we will denote the application of the operator 
d z ), W is a complex-valued function and W = C[W] is its complex conjugate 
function. As was shown in [TJ], [15] (see also [16], [17], [18]) equation fJTTT) is 
closely related to the second-order equation of the form 

(-A + i/)u = in Q (12) 

where v = (A/) // and u are real- valued functions. In particular the follow- 
ing statements are valid. 

Theorem 1 f73) / Let W be a solution of 07]) . Then u = Re W is a solution 
of ( [HP and f = Im is a solution of the equation 

(-A + ri)v = (13) 

where rj — 2 (^f^) 



2 
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Theorem 2 fI3) / Let u be a solution of (QI|) m a simply connected domain 
Q. Then the function 



v e ker A + v — 2 
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such that W = u + iv be a solution of < f77]j . is constructed according to the 
formula 

v = r^iifur 1 ^). (14) 

It is unique up to an additive term cf~ l where c is an arbitrary real constant. 

Given v G ker ^A + v — 2 (~p^ ^ 5 ^ e corresponding u G ker (A — u) 
can be constructed as follows 

u = -fA{if- 2 d^fv)) (15) 

up to an additive term cf. 

Thus, the relation between (JTT1) and (JT2l) is very similar to that between 
the Cauchy-Riemann system and the Laplace equation. Moreover, choosing 
v = and / = 1 we arrive at the well known formulas from classical complex 
analysis. Note that the potential rj has the form of a potential obtained after 
the Darboux transformation but already in two dimensions. 

For a Vekua equation of the form 

Wg=aW + bW 

where a and b are arbitrary complex-valued functions from an appropriate 
function space |28j, a well developed theory of Taylor and Laurent series in 
formal powers has been created (see [I], [5]). We recall that a formal power 
(a, zq, z) corresponding to a coefficient a and a centre Zq is a solution of 
the Vekua equation satisfying the asymptotic formula 

ZW(a,z ;z) ~a(z-z ) n , z - z . (16) 

For a rigorous definition we refer to [1], [5]. The theory of Taylor and Laurent 
series in formal powers among other results contains the expansion and the 
Runge theorems as well as more precise convergence results [2T] and a recently 
obtained simple algorithm [T7] for explicit construction of formal powers for 
the Vekua equation of the form ffTTT) . In section H] we show how this theory 
is applied for generalizing the second Euler theorem. For this we need the 
expansion theorem from [1]. For the Vekua equation of the form ffTTj) this 
expansion theorem reads as follows (for more details we refer the reader to 

mi)- 
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Theorem 3 Let W be a regular solution of < f77]j defined for \z — z \ < R. 
Then it admits a unique expansion of the form W(z) = Y^=o ^'"Vm z o', z ) 
which converges normally for \z — zq\ < R. 

Another result which will be used in the present work (section Hj) is a 
Cauchy-type integral theorem for the stationary Schrodinger equation. It 
was obtained in [T3] with the aid of the pseudoanalytic function theory. 

Theorem 4 (Cauchy's integral theorem for the Schrodinger equation) Let f 
be a nonvanishing solution of fiTty) in a domain Q and u be another arbitrary 
solution of (TB) in Q. Then for every closed curve T situated in a simply 
connected subdomain of Q, 

Re / d,£)dz + ilm [ f 2 d z (^)dz = 0. (17) 
Jr J Jr J 

3 The two-dimensional stationary Schrodinger 
equation and the complex Riccati equation 

Consider the complex differential Riccati equation 

^Q + \Q\ 2 = \ (18) 

where for convenience the factor 1/4 was included. We recall that v is a 
real-valued function. Together with this equation we consider the stationary 
Schrodinger equation f|T2|) . 

(-A + i/)w = (19) 
where u is real- valued. Both equations are studied in a domain Q C C. 
Theorem 5 Let u be a solution of [Ttfy . Then its logarithmic derivative 

Q = - (20) 
u 

is a solution of (Tig)). 

Proof. It is only necessary to substitute (T2U]) into (ITS]) . ■ 
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Remark 6 Any solution of equation l[T8\) fulfils / fiOj) . Indeed, the imaginary 
part of l[T8\) reads as follows 

d y Qi + d x Q 2 = 0. 
Theorem 7 Let Q be a solution of / flgj) . Then the function 

u = e A[Q] (21) 

is a solution of / figj) . 

Proof. Equation ( fT9|) can be written in the form 

(4c\d z - v)u = 0. 

Taking u in the form (12T|) where Q is a solution of (TTHT) and using the obser- 
vation that 

&(A[Q]) = &(A[Q]) = Q 

we have 

d# K u = d, (Qe A M) = e A ™ (ckQ + \Q\ 2 ) = \u. 

■ 

Observe that this theorem means that if Q is a solution of ([IB"]) then there 
exists a solution u of (fl9l) such that (!20|) is valid. Theorems [5] and [7| are direct 
generalizations of the corresponding facts from the one-dimensional theory. 

The following statement is a generalization of the one-dimensional fac- 
torization (BJ. 

Theorem 8 Given a complex function Q, for any real valued twice contin- 
uously differentiable function the following equality is valid 

~(A-v) ( p=(d s + QC)(d g -QC)(p (22) 
= (d s + QC)(d s -QC)<p 
if and only if Q is a solution of the Riccati equation fi~8\) . 
Proof. Consider 

{d- z + QC)(d z - QC)ip = l -Ay -\Q\\- 

from which it is seen that fl22l) is valid iff Q is a solution of (|18p . The second 
equality in fl22l) is obtained by applying C to both sides of the first equality. 
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4 Generalizations of classical theorems 



In this section we give generalizations of both Euler's theorems for the Riccati 
equation as well as of Picard's theorem. 

Theorem 9 (First Euler's theorem) Let Q be a bounded particular solution 
of (ESP- Then [W) reduces to the following first order (real-linear) equation 



Wz = Q W (23) 
in the following sense. Any solution of [W) has the form 

d z ReW 



Q 



ReW 



and vice versa, any solution of El3) can be ex-pressed via a corresponding 
solution Q of ffify as follows 

W = e A ^ + ie-^A [ie 2A ^c\e A ^ Q ^] . (24) 

Proof. Let Qo be a bounded solution of (|18p . Then by theorem [7] we have 
that there exists a nonvanishing real valued solution / of (fT9l) such that 
Qo = /*//• Hence ( l23l) has the form ffTT]) . Now, let Q be any solution of 
( TT8l) . Then again Q = u z /u where u is a solution of ( Tl9l) . According to 
theorem [2} u is a real part of a solution W of (ITTj) . Thus we have proved the 
first part of the theorem. 

Let W = u + iv be any solution of (123"]) (m = Re W). Then w is a solution 
of (fl9l) and by theorem [5] it can be represented in the form u = e A ^ where 
Q is a solution of (fl8|) . Then by theorem [2] (formula (fT4^) ) , W has the form 



Thus the Riccati equation ( 1181) is equivalent to a main Vekua equation of 
the form (TTTT) . 

In what follows we suppose that in the domain of interest Q there exists 
a bounded solution of (TT8l). 



Theorem 10 (Second Euler's theorem) Any solution Q of equation Iffi^) de- 
fined for \z — z \ < R can be represented in the form 

V°° O e A ^ 

Z— /n=0 
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where {Q n }^Lo ^ s ^ e se ^ of particular solutions of the Riccati equation l[T8\) 
obtained as follows 

d z ReZ^(a nj z ,z) 
W) ~ ReZ(n)(a n ,z ,z) ' 

Z^ n \a n , Zo, z) are formal powers corresponding to equation [2B) and both se- 
ries in If25*\) converge normally for \z — zq\ < R. 

Proof. By the first Euler theorem we have 

d?ReW 



Q 



ReW 

where W is a solution of ( 123|) . From theorem [3] we obtain 

1 } EZo^Z {n K*n,z -z) ■ 
Every formal power Z^ n \a n , z ; z) corresponds to a solution of ([TBI) : 

d z ReZ^{a m z ;z) 



Qn(z) 



ReZ( n )(a n ,z ;z) 



or Re Z^(a n , z ; z) = e A ^ Qn ^ from where we obtain (1251) . ■ 

In the next statement we give a generalization of Picard's theorem in the 
form (J6]). 

Theorem 11 (Picard's theorem) Let Qk, k = 1,2,3,4 be four solutions of 
(ED). Then 

dziQt - Q 2 ) + 1% Im^Qa) c\(Q 3 - Q 4 ) + 2t Im(Q 3 Q 4 ) 



Qi — Q2 Q3 — Qa 

dz(Qi - Qa) + 22 ImiQM c\(Q 3 - Q 2 ) + 2i lm(Q 3 Q 2 



Qi — Qa Q3 ~ Q2 

Proof. Obviously, 

(Qi + Qa) + (Q 3 + Q2) - (Qi + Q2) ~ (Qs + Qa) = 0. 



0. 
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Multiplying each parenthesis by 1 = (Qi — Qj)/(Qi — Qj) we obtain the 
equality 

(Q1 + Q4XQ1-Q4) , (g 3 + Q 2 )(Q3-Q 2 ) 



(Qi - Qi) (Q 3 - Qi) 

(Qi + Q 2 ) (Qi - Q2) (Q 3 + Qd (Qs - Q±) 



(Qi - Q2) (Q 3 - Qa) 

Using 

(Qi + Q 3 ){Qi - Qj) = dz{Qj - Qi) - QiQj + Q& 

the result is obtained. ■ 



0. 



5 Cauchy's integral theorem 

Theorem 12 (Cauchy's integral theorem for the complex Riccati equation) 
Let Qo and Qi be bounded solutions of [Tfy) in Q. Then for every closed curve 
T lying in a simply connected subdomain of Q, 

Re J {Q x - Qo) e A[Ql ~ Qo] dz + i Im J {Q x - Q ) e A[Ql+Qo] dz = 0. (26) 

Proof. From theorem [7] we have that / = e A '^ and u = e A ^^ are solutions 
of (Tl9l) . Now, applying theorem H] we obtain 

Re / d t £)dz + i1m [ f 2 d z (^)dz = 
Jr f Jr J 

for every closed curve T situated in a simply connected subdomain of Q, 
which gives us the equality 

Re J d z (e A[Ql - Qo] )dz + ilm J e 2A[Qo] d z (e A[Ql - Qo] )dz = 0. 

From this we obtain the result. ■ 

As a particular case let us analyze the Riccati equation f|T8l) with v = 
which is related to the Laplace equation. If in fl26l) we assume that Q = 0, 
then (1261) takes the following form 

J^Q ie A[Q ' ] dz = 0. 
11 



This is obviously valid because if Qi is another bounded solution of the 
Riccati equation with v = 0, then according to [7] we have that u = is 
a harmonic function and the last formula turns into the equality 

u z dz = (27) 

jt 

(u z is analytic). 

Now, if in (12T)1) we assume that Qi = 0, then (12BT) takes the form 
Re J Q e- A[Qo] dz + ilm J Q e A[Qo] dz = 0. 

Rewriting this equality in terms of the harmonic function / = e A ^ ^ we 
obtain 

fzdz 
f 2 

which taking into account (I27|) becomes the equality 



Re / ^V + zIm / fzdz = 



or which is the same, 

Re J d z (jj dz = 0. (28) 

This equality is a simple corollary of a complex version of the Green- Gauss 
theorem (see, e.g., [271 sect. 3.2]) according to which we have 

d z ( \ J dz = djd z ( - J dxdy. 



For / real the right-hand side is real- valued and we obtain ( 1281) . 



6 Conclusions 

We have shown that the stationary Schrodinger equation in a two-dimensional 
case is related to a complex differential Riccati equation which possesses 
many interesting properties similar to its one-dimensional prototype. Be- 
sides the generalizations of the famous Euler theorems we have obtained 
the generalization of Picard's theorem and the Cauchy integral theorem for 
solutions of the complex Riccati equation. The theory of pseudoanalytic 
functions has been intensively used. 
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